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£NJ , Abstract 

The semiclassical limit of a 4-simplex amplitude for a spin foam quantum gravity model 
with an Immirzi parameter is studied. If the boundary state represents a non-degenerate 
4-simplex geometry, the asymptotic formula contains the Regge action for general relativity. 
A canonical choice of phase for the boundary state is introduced and is shown to be necessary 
to obtain the results. 
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1 Introduction 

u . 

A key step in understanding the semiclassical regime of a spin foam model [U El E] in dimen- 
sion n is the analysis of the asymptotic behaviour of the n-simplex amplitude that defines the 
model. In fact, the discovery by Ponzano and Regge that the 6j symbol of recoupling theory 
contains the Regge action in its asymptotic behaviour established it as a model for 3D quantum 
C*"*- . gravity [3]. Similar asymptotic analysis of the 4-dimensional models [5] was initially performed 

by Barrett, Williams and Steele [fU [7], and formed the basis of investigations of the graviton 
propagator structure of these models [HJ EJ [10] . This latter analysis showed a definite incom- 
patibility between the lOj symbol and a boundary structure given by loop quantum gravity-like 
geometry. Consequently a host of new 4-dimensional models were developed. The first such 
models, were due to Engle, Pereira and Rovelli [UJ. Meanwhile Livine and Speziale introduced 
coherent states to the analysis and definition of spin foam models |12j . and suggested a way to 
construct new models in [13]. In parallel, Preidel and Krasnov defined and developed a model 
along these lines in p3j. A refined version of the original EPR model was published in [15] 
together with Livine. The FK model and the EPRL model depend on the Immirzi parameter 
7 and are identical for 7 < 1. An initial exploration of the asymptotics of the FK model for 
manifolds without boundary was undertaken in [16j . 

Our thanks are due to Carlo Rovelli, who encouraged us to study this problem, and for 
hosting JWB and WF in Marseille, where the methods and most of the results were presented 
in detail on 19th November 2008 in a five-hour seminar. The result presented in Marseille was a 
complete derivation of the term in the asymptotic formula which gives the cosine of the Regge 
action for 7 < 1, and is described by Alesci, Bianchi and Rovelli in [T7], and also in the talk by 
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Alesci at PI |18j . The derivation of the phase part of this term is also confirmed by the work of 
Freidel and Conrady in [19j . 

A spin foam model is a procedure to compute an amplitude from a triangulated manifold 
T with n-simplices A n coloured by representation theory data. In four-dimensions, such an 
amplitude is typically of the form 

z ^ = e n n o n f*o» o- « 

hp A 2 A 3 A 4 

where / n are weights assigned to the n-simplices of the triangulated manifold, and p and i 
respectively denote the assignments of unitary, irreducible representations to the 2-simplices, 
and intertwining operators to the 3-simplices of T. The model is specified by the choice of 
representation assignments, the vector space of intertwining operators l, and weights f n . 

The Euclidean EPRL is a specific spin foam model [15] built out of the representation theory 
of the spin group G = Spin(4). In this paper we analyse in detail the large spin structure of the 
four-simplex amplitude of the model. We begin by defining the EPRL model in section [2] In 
section [3] we discuss the general well-posedness of the analysis to be performed here and discuss 
the geometry of the three dimensional boundary determined by the chosen intertwiners l. In 
section 0] we state the main result of the full vertex amplitude, which we then prove in section 
[5] We conclude by discussing some further aspects of the results obtained in section [6] 

2 The EPRL four-simplex amplitude 

The input data for the 4-simplex amplitude is a spin k £ {0, ^, 1, . . .} for each triangle of the 
4-simplex and an SU(2) intertwiner t for each tetrahedron. From t, a Spin(4) intertwiner i is 
constructed, and then these Spin(4) intertwiners are glued together in the standard fashion to 
construct an amplitude (a complex number) for this data. The only other input required is the 
Immirzi parameter 7, which is a constant. 

Firstly, a precise definition of t is required. For a given tetrahedron, one has to choose an 
ordering of its four faces, e.g., by numbering them with 1,2,3,4. Then the SU(2) intertwiner 1 
is an element of 

4 

Hom su(2) (C,(g)y fc J, 
i=l 

where the spaces are tensored together in the order <8> Vfc 2 <8> Vk 3 <8> V]~ A ■ This ordering convention 
is used throughout. 

Of course the spaces constructed using different orderings are easily related by an action of 
the permutation group. We use the binor category of representations [201 [2T| [22] throughout 
the paper. In this category the crossing diagram is fermionic, which means that the crossing 
of two lines of odd spin gives a factor of —1. For example, the map T4 1 (g) Vk 2 — > Vk 2 <g> Vk ± is 
x (g> y — ► (— l) 4fclfc 2y x _ Spin network diagrams in this category can be evaluated using the 
Kauffman bracket [23] specialised to Kauffman's parameter A = — 1. The binor calculus has the 
convenient feature that the framing of a curve does not affect the evaluation. 

The i are constructed as follows. Let (tt&, Vk) and (7Ty- Vu- respectively denote the 
unitary, irreducible representations of SU(2) and Spin(4) = SU(2)_ x SU(2) + . There exists an 
injection 

4 4 
<P : Hom su(2) (C, (g) V h ) Hom Spin(4) (C, (g) V^-^) 

i=l i=l 

I 1 ► <j)(i) := l, (2) 
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■fclfc2fc 3 fc4 



Figure 1: The Spin(4) intertwiner i. 



embedding the vector space of SU(2) intertwiners into the vector space of Spin(4) intertwiners. 

Explicitly, 4> is constructed by using the Clebsch-Gordan interwining maps 5 ' : — ► 
Vj- ® Vj+ injecting the SU(2) representation V*. into the highest (resp. lowest) diagonal SU(2) 
subgroup factor k = j + + j~ (resp. k = j + — j~) of Vu- j+\ = Vj- £ED in the 7 < 1 (resp. 
7 > 1) case. The labels and k are related via the Immirzi parameter by 

^ = ^±7^. (3) 

These relations of course constrain the values of k so that the are always half integers; 
specifically if 7 = p/q is written in lowest terms, then k has to be a multiple of either q/2, or q 
in some cases. 

The Spin(4) intertwiner 1 is then obtained as follows 

r 4 .--+ 

1 := 0(0 = / dG fa- ® 7T.+ )(G) o (g) Jl o ^ fe ^3*4 > (4) 

^Spin(4) ^ Ji 111 1 

where the notation G = (X~ , X + ) is used (see figure [TJ. The group integration ensures that the 
resulting object is Spin(4)-invariant, i.e., is an element of Homg pin (4)(C, (S>t=i Vu~ 

The dynamics of the EPRL model is encoded in the four-simplex, or vertex amplitude 
constructed by contracting the specified intertwining operators associated to each of the five 
tetrahedra A3 of the four-simplex A4. Labelling these tetrahedra by a = 1, 5, the ten triangles 
A2 of A4 are then indexed by the pair ab of tetrahedra which intersect on the triangle. There 
are two SU(2) group elements (X~, X£) and one SU(2) intertwiner t a for each tetrahedron, and 
three (7-related) SU(2) representations k a b and (j^,^) for each triangle. The intertwiner t a 
lies in the space 

Hom su(2) (C, (g) V k J, 

b:b^a 

with the tensor product ordered by the numerical order of b, and k a b = k^a- 

The amplitude G C is defined by forming a closed spin network diagram from this data. 
The five intertwiners (vertices) i a are tensored together and then the free ends are joined pairwise 
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according to the combinatorics, i.e. the edge a of vertex b is joined to edge b of vertex o. This 
is done using the standard 'e inner product' of irreducible representations of SU(2), denoted 



e k : V k ® V k -> C. 

This is defined by a choice of the two-dimensional antisymmetric tensor e for SU(2) spin 1/2, 
and extended to arbitrary spin by tensor products of e. This inner product is represented in the 
spin network diagram as a semicircular arc. To combine the Spin(4) intertwiners L a , one regards 
each vertex as an SU(2) spin network and uses one e inner product to connect the j + edges and 
a second e inner product to connect the j~ edges. The resulting closed diagram is evaluated 
using the binor calculus rules for all crossings. (Note: one could also use the rule that there is 
no sign for a crossing of a + line with a — . This makes at most a difference of an overall sign to 

AO 

This yields a formula 



h = (-i) x / X[dxZdx- ( 0^ ) o ((g) i a 



(5) 



where the propag ators K ab : V k * b ® V k ^ b -> C are defined by 



and (— l) x is the sign defined by the diagrammatic calculus of spin networks. Notice that as is 
linear in the t we can use unnormalized intertwiners and push the normalization into /a(£) 2 = — , 
the asymptotic behaviour of which can then be analysed independently. 

The aim of this paper is to study the large spin behaviour k a b -> oo of the four-simplex 
amplitude f&. Our strategy is to approximate the integral formula of the amplitude using 
(extended) stationary phase. At this level, the problem is not well posed because the scaling of 
the SU(2) intertwiners is not yet defined. Solving this problem naturally leads to a reformulation 
of the integral formula to an exponential form which is particularly well suited to asymptotics. 

As a final remark, although we will not use this in the following, using the SU(2)-invariance 
of the Clebsch-Gordan maps, one can set one of the two group arguments of the propagator, 
say the left handed part, to the identity. In this case, the amplitude © becomes the Feynman 
evaluation associated to a tensor field theory over S 3 = SU(2). The 'matter fields' are identified 
as sections of the vector bundle E = P x k V k associated to the trivial principal bundle P = 
Spin(4) = S 3 x SU(2) with base manifold 5 3 and structure group SU(2) via the representation k. 
The internal indices of the propagators of the matter fields are contracted at the vertices of the 
diagram using the SU(2) intertwiners t, and the amplitude with X~ = 1, is the Feynman 
evaluation associated to the complete graph with five vertices. 

3 Posing the asymptotic problem 

The asymptotic large spin limit can be investigated by scaling each spin simultaneously, replacing 

then fixing numbers k ao , and taking the limit A — > oo through the natural numbers. However 
to specify the boundary data correctly for each A, one needs to 'scale' the corresponding inter- 
twiners. To do this, the intertwiners are specified in a particular way which gives a sequence 
of intertwiners in the corresponding spaces. The intertwiners are specified with coherent states, 
following the construction in 
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3.1 Coherent states 



Let n £ S 2 be a unit 3-vector. Then a coherent state a S V& in direction n is a unit vector 
satisfying 

(L.n) a = ika, 

with L the standard anti-hermitian rotation generators in the k representation and the dot '.' 
denoting the 3d (Euclidean) scalar product. The coherent state has the maximal spin projection 
along the n axis. 

One way of thinking of coherent states is to pick an arbitrary coherent state T(n) for each 
unit vector n. Then any other coherent state is a phase factor e ie times one of these standard 
ones, a = e td T(ii). This information is displayed in bra-ket notation as 

a = \k, n, 9). 

In fact the choice of section V does not play any role in the following. 

The notation is shortened in the following ways. If the parameter k is omitted, then the 
fundamental representation is used; thus |i,n, 0) = |n, 0). Further, the parameter 6 will be 
omitted later when the phase is obvious from the context. 

Embedding Vk in <S> 2k Vi/2 shows that 

® 2k \n,9) 

is a coherent state for n and Vk ■ Therefore the Hermitian inner product of coherent states obeys 

\(k,n 1 ,e 1 \k,n 2 ,0 2 )\ 2 = U(l + ni.n 2 )J , (6) 
where | . | 2 denotes modulus square. 



3.2 Coherent tetrahedra 

To construct a coherent intertwiner for a given tetrahedron, the idea is to associate a coherent 
state to each one of its triangles. The geometrical picture is that the coherent state \k, n, 6) then 
carries the interpretation of the normal of length k and direction n to the associated triangle, 
plus a phase factor. 

Furthermore, we want to describe tetrahedra with three-dimensional rotational symmetry 
and the coherent intertwiners are thus constructed by integrating over all spatial directions the 
tensor product of four coherent states 

r 4 

^ lfc2fc3fc4 (ni,n 2 ,n 3 ,n 4 ) = / dh (g) h\k h n h 0$. (7) 

^SU(2) i=1 

These intertwiners were introduced by Livine and Speziale [12J, who gave an asymptotic formula 
for their normalisation. 

According to the 'quantization commutes with reduction' theorem of Guillemin and Sternberg 
[25], the space of intertwiners is spanned by the t determined by vectors satisfying the closure 
constraint 

%ni + k 2 n 2 + /c 3 n 3 + /c 4 n 4 = (8) 

Therefore in this paper the coherent tetrahedron states are always taken to satisfy this condition. 
The condition also implies there is a tetrahedron t € M 3 in Euclidean space, with the standard 
metric, which has these four vectors as the outward face normals and triangle areas equal to 
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ki |26j . Throughout this paper attention is restricted to the cases where the tetrahedra are 
non-degenerate. This means that the analysis excludes a few extremal cases. 

The tetrahedron t is uniquely specified by the four vectors ni , 112 , 113 , 114 , up to parallel trans- 
lation. Therefore the tetrahedron has a definite orientation. The coherent state i is averaged by 
the action of SU(2) on the coherent states, which amounts to an action of SO(3) on the tetra- 
hedron t, i.e. rigid rotations which preserve the orientation. The geometry of the tetrahedron 
which is invariant under these rotations is a metric and an orientation. The coherent state is a 
quantum version of the geometry of this tetrahedron. 

3.3 Coherent boundary 3D manifolds 

3.3.1 Boundary data 

For the asymptotic formula in this paper, the boundary manifold is just the boundary of the 
4-simplex, i.e. a three-sphere triangulated with five tetrahedra. However the considerations of 
this section apply to more general cases. So consider E to be a triangulated closed 3-manifold. 
To specify a coherent state for the whole of E, the information required is: i) a choice of spin 
k for each triangle, ii) for each tetrahedron r, a choice of coherent state \k,n,9) for each of its 
four triangles. Clearly the k used in ii) are the ones picked in i). Such a choice B = {k,n} for 
the whole of E is called boundary data for the manifold S. 
Then the boundary state specified by this data is 

1>(k,n)= (g) ^^(m.n^na.iu). 

tetrahedra 

However the phase of this state vector is arbitrary, as the phase of each coherent state has not 
been specified yet. In the following, it is shown that for an important subset of boundary data, 
there is a canonical choice of phase. 

3.3.2 Regge-like boundary data 

For each tetrahedron r a C E let <p a : r a — > M 3 be the affine linear map such that (p a {T a ) = t a is 
the geometric tetrahedron defined above, i.e., such that is the outward normal of t a in the 
direction of the neighbouring tetrahedron b. 
Boundary data will be called Regge-like if 

• There is a metric g-% on E (as in Regge calculus) such that the metric restricted to the 
tetrahedron r a is the pull-back of the standard metric 5 on M 3 with <j> a , i.e. gT,\r a = <?^- 

• There is an orientation Or of E such that the orientation Or Ta of each r a is the pull-back 
of the standard orientation on M 3 with (j> a . 

Of course these conditions can only be satisfied if E is an orientable manifold. 

3.3.3 Regge states 

For Regge-like boundary data there is a canonical choice of phase for the state, which will now 
be described. 

The construction needs the standard antilinear structure map for representations of SU(2), 
J : T4 — > Vk ■ This is defined by 

e^(a, a ) = {J a\a ), 
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the left-hand side being the epsilon-inner product and the right hand side the Hermitian inner 
product. It obeys Jg = gJ for all g € SU(2), J 2 = (-l) 2k and (Ja\Ja') = (a\a'). Since 

J(zn • L) = —(in • L) J, 

the map J takes a coherent state for n to a coherent state for — n. 

Recall that the SU(2) coherent state \k,n,9) represents a vector n of fixed length in IR 3 , 
interpreted as describing the normal vector to a triangle lying in the plane n- 1 orthogonal to n 
in IR 3 . This triangle is only defined up to U(l) rotations in the corresponding plane (this is the 
phase information), and the coherent state \k, n, 9 + <p) describes the same triangle but rotated 
by an angle 2tp in the n -1 - plane. 

Let T a and t& denote two tetrahedra in S which share a common triangle, A a b = r a n r&. The 
important point is that for Regge-like boundary data the geometries of the common triangle 
agree. This means that there is a unique element g a b 6 0(3) which, together with a translation, 
maps the image of this triangle (j) a {A a b) m ta congruently to the corresponding triangle fa(A a b) 
in tf, and takes one outward-pointing normal to minus the other one, 

g a b O fa{A ab ) = fa(A ab ) 

Qab^-ab = -n 6a . (9) 

In fact, g a b is in S0(3), due to the following argument. Consider any two linearly independent 
vectors v, w in A ab . The element g a b is defined to be the linear transformation that maps </> a ( v ) 
to fa(v), fa(w) to fa(w) and n a b to — nj, a . However, due to the fact that (f) a and fa are 
orientation-preserving, these vectors form frames that have the same orientation in IR 3 . Thus 

gab g S0(3). 

These group elements can actually be lifted to g a b G SU(2), using a choice of spin structure 
for the manifold S. Given a spin structure, the Levi-Civita connection for the metric lifts to 
a spin connection on the spin bundle. As the manifold is not smooth everywhere this is a slight 
extension of the usual notion, but is determined in a fairly straightforward way by smoothing 
the conical singularities on the 1-skeleton of S (the edges and vertices) and applying the usual 
definition to the smoothing. In fact the holonomy of a loop around a conical singularity is not 
just a rotation matrix, but in fact is a rotation angle (which deforms continuously to zero when 
the loop is shrunk to a point in a smoothing of the cone) . In particular this determines uniquely 
the holonomy in the spin group around a conical singularity. 

The Levi-Civita connection determines a parallel transport operator from the tangent space 
of one tetrahedron r a to another, r& 

u a b ■ Tr a — > Tn- 

The link with the matrix g a b is to regard it as the parallel transport operator in bases provided 
by the tangent space frames (f> a and 4>b f° r each tetrahedron, i.e., 

9ab = faUablfia 1 , 

disregarding the translation part of these maps. Concretely, this can be regarded as the parallel 
transport along a dual edge, with the frames 4> a and fa at the corresponding dual vertices. 

Given a choice of spin structure for S, there are two spin frames for each tetrahedron that 
cover the given tangent space frames. Choosing one of these for each tetrahedron then defines 
the g a b £ SU(2) as the parallel transport operators for the spin connection using these spin 
frames as bases. 
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Figure 2: The propagator V a b for a single edge. 



The choice of phase for the boundary state V is given by picking the phase of a. a b = 
\k a b, n ab,8ab) for T~a to be arbitrary, and then fixing the phase of the state for the corresponding 
triangle in r& to be 

®ba = 9abJot a b- 

In the definition of ip, both of these states appear; ip is an integral over a product of terms 

Oiab ® gabJoiab-, 

one for each triangle. Since J is antilinear, a change of phase of a cancels between these two 
factors of the tensor product. This provides a canonical choice of the phase for the boundary 
state t/j. Regge-like boundary data together with this choice of phase of the state ip is called a 
Regge state. 

Although the g a \, £ SU(2) depend on the choice of a spin frame for each tetrahedron, changing 
the choice of spin frame does not affect the Regge state. This is because changing the spin frame 
for tetrahedron a to the other possible choice results in the simultaneous change g a b — > — g a b for 
all tetrahedra b neighbouring a. Since the sum of the spins on the face of tetrahedron a is an 
integer, the action of ®b{~9ab) on the tensor product of the coherent states is the same as the 
action of <g>b9ab- Thus the Regge state is the same for either of the two possible choices of spin 
frame at any given tetrahedron. It depends only on the chosen spin structure. Of course, for 
the case of a 4-simplex, the boundary X = S 3 has only one spin structure and this dependence 
is of no consequence. 

It is worth noting that, in this context of general boundaries X, Regge states will behave well 
under gluing. In other words, gluing simplexes together with the canonical choice of phase will 
produce a formula for the partition function for a general manifold with the canonical choice of 
phase on the boundary. 

From now on the notation for coherent state vectors will be shortened to \k, n), omitting the 
phase parameter. The notation \k,— n) means J\k,n). 
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3.4 Exponentiated expression for the four-simplex amplitude 

The general considerations of the previous sections are now applied, for the rest of the paper, to 
the case of a 4-simplex a. The boundary data is specified on the simplicial 3-manifold £ = do. 
Using the coherent states framework, the four-simplex amplitude is 



h = (-i) x ' / TT dG a / TT dh* TT r ab , 

7 S pin(4)5 Y ^SU(2)5 Y Yb 



(10) 



where the coherent propagator Pat is now given by 



ab 



^abi 



lablab 



-n a »kU^)^+ « r jX*) , 3+ jX+) C kat 



Kk ab (h b )\k ba ,n ba ), 



(11) 



using the notation X^ b := (X^)~ 1 X b iz , and C k -^- + : \% — ► V}_ (8) Vj + the reflected spin network, 
as shown in figure [21 This is proved starting from ([5]), rotating some of the Clebsch-Gordan 
coefficients using the epsilon inner products, then converting these inner products to the Her- 
mitian inner product, and also flipping the order of + and — lines, obtaining further factors of 
— 1 which are absorbed into the definition of 

The next step is to obtain an exponentiated version of the amplitude as a means to enter 
the framework of (extended) stationary phase. This is realized through a reformulation of the 
propagators. 

The first remark in order is that the integration over SU(2) in equation (1101) at each vertex 
can be absorbed into the Spin(4) integration because of the invariance of the Clebsh-Gordan 
maps. Then, the idea is to use the exponentiating property of the coherent states 



I A;, n) 



In) 



(g)2fc 



(12) 



to reduce the propagator to a product of two propagators in the fundamental representation of 
SU(2) to the power 2j± respectively. 

For the next step, we need to treat the 7 < 1 and 7 > 1 cases separately. 



7 < 1 case. For 7 < 1, the Clebsch-Gordan coefficient C 1 k ah ^ ab injects into the highest spin 
subspace k ab = + j~ b . Considered as a spin network, is easy to see that the symmetrizers 
on the and j~ b edges can be absorbed into the symmetrizer on the k ab edge because of the 
stacking property of symmetrizers, see figure O The remaining symmetrizer now acts on the 



jab + 1 




ab 



Figure 3: The three-valent intertwiner cl ab:lab for the case 7 < 1 showing how the projection to 

K ab 

the highest spin subspace k ab makes two of the symmetrizers redundant. 
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coherent states \k ab ,n ab ) but since these are defined as the symmetrized tensor product of spin 
half coherent states \n ab ) this final symmetrization can also be ignored. We can now use the 
exponentiating property of the coherent states to further split the propagator into a product of 
terms in the fundamental representation. We obtain the following expression for the propagator 



P. 



7<1 



ab 



n ab \X , \n ba ) 2j ob (- na6 |Xi|n 6a ) 



■' a 



The four-simplex amplitude can thus be re-expressed as 

/4=(-i) x ' / n^ e57<i > 

JSpin(4)5 a 

with the action given by 

S 7<i = Yl 2j ab ln(-n a6 |X- 6 |n 6a ) + 2j+ In (-n ab \X+ |n 6a ). 



(13) 



(14) 



(15) 



a<b 



Notice that this action is in general complex. The logarithm of a complex number is only defined 
up to a multiple of 27ri, we can safely neglect this factor as it will not affect the stationary points 
and when it appears in the action it is exponentiated. 

7 > 1 case. The 7 > 1 case works analogously but is more complicated as the Clebsch-Gordan 
coefficient C J k abJab now injects into the lowest state k ab = — j~ b , see figured! 




«o6 



•Jab Jab 




3ab 



Figure 4: The three-valent intertwiner C J k ab J ab for the case 7 > 1 projecting to the lowest weight 
state k ab =j+-j+. 



The symmetrizers on the k ab and j~ b edges can be absorbed into the symmetrizer on the ft 
edge. We remove the remaining symmetrizer by using the resolution of the identity in terms of 
SU(2) coherent states in the j^ b representation 



lj = dj / dn\j,n){j,n\ (16) 

JSU(2)/U(1) 

where an arbitrary choice of phase is made for each coherent state, dj = 2j + 1, and the integral 
measure is normalised to 1. 
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With these insertions, all of the symmetrizers can be absorbed into coherent states and the 
propagator can be divided into terms in the fundamental representation as before 

Kb' = (- 1 ) 2ia ~ 6d ?+/ dm ^ j dm ba (-n ab \m ab ) 2k ">(m ab \X+\rn ba ) 2 ti 

x(m ab \X- b \ui ba ) 2jab (m ba \n ba ) 2k ^. (17) 
The action is given by 

5 '7>i = 2 iab ln ( m ab\X+ b \m ba ) + 2j~ b In (m afe |X~ |m fea ) 

a<b 

+2A; a6 ln(-n a ; ) [m afe } + 2k ab ln (m 6a [n 6a ) (18) 

Note that we could just have easily decomposed the propagator using only the m ab , for a < b. 
However, using the second resolution of the identity gives a formula that treats the ± sectors in 
a more symmetric way. 

The j~ term is easiest to understand in terms of the graphic notation. We first flip over 
the left bend in the diagram. This causes the factor of (— l) 2 - 7 in equation (|17p . Then we can 
straighten out the "5"' by using standard graphical calculus, which leads to a transposition of 
the group element. Written as equations in the fundamental representation this reads €jiX 3 k e kl = 
—€ijXle kl = —(X -1 )', using the convention ejj-e jfe = 5f. In terms of the standard inner product 
this part of the diagram then leads to the matrix element 

-(mhaKX^^lmah) = -(m ab \X~ b \m ba ). 
Thus the amplitude for 7 > 1 is given by 

U = ("l) x ' TT(— l) 3 ^*-* d\ [ T\dG c [ dm ab dm ba e s T>i. (19) 
3.4.1 Symmetries of the action. 

The actions (I15D and (I18p admit two types of symmetry which will be important later. 

• Continuous. A global Spin(4) transformation characterised by an element (Y~,Y + ) in 
Spin(4) acting on each X^ by X^ — > Y^X^. This determines a rigid motion of the whole 
4-simplex. 

• Discrete. At each vertex a the transformation X+ — ► — X+ leaves a factor (— 1)^,M<» 2j °f'. 
Since the j^ b sum to an integer at each tetrahedron a, this factor will equal one. This 
symmetry is a consequence of the fact that SU(2) is the double cover of SO(3). There is 
a similar symmetry for the minus sector, X~ — >■ — X~ . 

4 Asymptotic formula 

To study the semi-classical approximation of the theory we consider the large spin limit, and 
study the asymptotics. We start by scaling all ten spins by a constant parameter k ab — ► Xk ab . 
Then the scaled boundary state used is 

?P\ = 1p{^ab,Xk a b)- 

If the boundary data is Regge-like then the Regge state is used here. Otherwise the phase of 
the state is undetermined. 
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The asymptotic behaviour of the 4-simplex amplitude f± depends on the geometric character 
of the boundary data. Recall that it is assumed throughout that for each tetrahedron r a , the 
four vectors n a &, b ^ a, are the outward normals to the faces of a non-degenerate tetrahedron 
t a in three-dimensional Euclidean space. 

There are two important types of boundary data for a 4-simplex. The first is for Regge-like 
boundary data where the boundary metric is the boundary of a non-degenerate four-dimensional 
Euclidean geometry for the 4-simplex. In other words, the t a geometries fit together in four- 
dimensional space. 

The second type which is important is called a vector geometry. The boundary data (not 
necessarily Regge-like) is said to be vector geometry if there is unit vector v a {, £ M 3 assigned to 
the ab-th triangle satisfying 

^ KbVab = 
b:b-^a 

V Q fe = -V ba 

and for each tetrahedron, the n a fe are congruent to the corresponding v ab (i.e. n a b = h a v a b for 
h a € SO(3)). 

Examples of vector geometries are determined by a non-degenerate 4-simplex. This gives 
two distinct vector geometries, given by the self-dual and anti-self-dual parts of the triangle 
bivectors. Another example of a vector geometry is given by a linear map of a 4-simplex into 
R 3 . 

In the case of a non-degenerate 4-simplex, the geometry of the simplex is uniquely defined 
by the geometry of the boundary, and hence by the boundary data. For this result in convex 
rigid geometry see for example Theorem 4.10 in [28]. Therefore there are well-defined dihedral 
angles < Q a b < it for each triangle. These are the angles between the outward normals to the 
tetrahedra embedded in R 4 . These angles are used in the asymptotic formula. 

A formula is asymptotic if the error term is bounded by a constant times one more power of 
A -1 than that stated in the asymptotic formula. 

Theorem 1. (Asymptotic formula) Given a set B = {n^, /c a b} a ^{, of boundary data, then in 
the limit A — > oo 

1. If B is a non- degenerate 4-simplex geometry, and ip\ is the associated Regge state, then 

2Nl_COS I A 7 J>a6@a6j 
V a<b J 

+ Nl + exp ( iX Y k a bQab J + Nl_ exp I -iX ^ k ab Q 



~ (-i)*' (^) 12 



>ab 

a<b / \ a<b 



(20) 



The numbers Nl_, Nl + , A2_ are independent of X and are defined below. 



2. For a set of boundary data which is not a non- degenerate J^-geometry but forms a vector 
geometry the asymptotic formula is 

The number N is independent of X and is defined below. 

3. For a set of boundary data that is neither a non- degenerate ^-geometry nor a vector geom- 
etry, the amplitude is supressed and does not contribute for large X. 

h(^x) = o(X- M ) VM (22) 
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5 Proof of the asymptotic formula 



In this section, we provide a proof of the above theorem. We first review the notion of (extended) 
stationary phase, showing how the asymptotic formula is dominated by the critical points. Next, 
we calculate these points and relate the non-degenerate critical points to geometric four-simplices 
in M 4 . Finally, we evaluate the integrand of the integral formula on the non-degenerate critical 
points and show how to recover terms related to the Regge action for the corresponding four- 
simplex. 

5.1 Extended stationary phase 

To study the semi-classical approximation of the theory, we are interested in the limit of ex- 
pression (fT4"j) or ([19]) when A — > oo. Our strategy is to use extended stationary phase methods, 
that is, stationary phase generalized to (non purely imaginary) complex functions. The general 
approximation scheme is stated as follows [24] . 

Let D be a closed manifold of dimension n, and let S and a be smooth, complex valued 
functions on D such that the real part HeS < 0. Consider the function 

/(A) = I dxa(x)e xs ( x \ (23) 
Jd 

The Hessian of S is the n x n matrix denoted H. The stationary points are assumed to be 
isolated and non-degenerate; det H ^ 0. 

In the extended stationary phase, the key role is played by critical points, that is, stationary 
points for which ReS = 0. If S has no critical points then for large parameter A the function / 
decreases faster that any power of A -1 . In other words, for all N > 1: 

/(A) = o(X- N ), (24) 

If there are critical points, then each critical point contributes to the asymptotics of / by a 
term of order \~ n l 2 . For large A the asymptotic expansion of the integral yields for each critical 
point 

t) jmm eXSi '° ){1+0{im - (25) 

At a critical point, the matrix —H has a positive definite real part, and the square root of 
the determinant of this matrix is the unique square root which is continuous on matrices with 
positive definite real part, and positive on real ones. 

If S admits several isolated critical points with non-degenerate Hessian, we obtain a sum of 
contributions of the form (|25|) from each of them. So, to compute the dominant terms in the 
asymptotics of (1191) for large spins, we need to find the stationary points of the action S and 
restrict to those with zero real part. 

5.1.1 Stationary points 

We separate the variational problems of the 7 < 1 and 7 > 1 sectors. 

7 < 1 case. The stationary points of 5' 7 <i, determined by varying the group variables G = 
(X~ , X + ), are the same as the stationary points of 

I 1<x •- e ^<i = J](_n ab |X a -|n fea ) 2 ^ (-n ab \X+ |n ba ) 2 ^>. (26) 

a<b 
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The variation of an arbitrary SU(2) group element X, and the associated variation of its inverse 
X^ 1 are given by 

5X = T o X 5X- 1 = -X- 1 o T, (27) 

where T = \iT r a r is an arbitrary element in the Lie algebra su(2), and a r are the Pauli matrices, 
with eigenvalues ±1, satisfying 

a r a s = 6 rs l + ie rst a t . 
The variation of I leads3 to two complex 3d vector equations for each tetrahedron a: 

5I y<1 = & Vo £ {1, 5}, Y, lab ^ab = °> (28) 

b: b^a 

with 

, _ {-n ab \(Xt)- l aX±\n ba ) 

ab (-n ab \(X±)-iX±\n ba ) ' i29j 

where we have introduced a bold letter notation for three-dimensional vectors: V = (Vi, V2, V3) 
and a = (oi, o"2, (J3). The minus signs arising from (I27p can be absorbed by noting that = 
— This can be proved by using the fact that the formula is the ratio of two epsilon inner 
products, and that whilst e(ga, (3) = e{a,g~ l (3) for a group element g, for a Lie algebra element 
T one has e(Ta,(3) = —e(a,T/3). 

As remarked above, each of these equations is complex. To extract a real and an imaginary 
part will need a bit more work. Firstly, we use the fact that the action of an element of SU(2) 
on a coherent state produces a new coherent state (temporarily putting the phase information 
back in the notation) 

\<b^ab) --=XtKb, 



with the three-vectors 



Xa^ab- 



In this formula the SU(2) element acts on a 3-vector via the homomorphism to SO(3). We will 
use this notation for the remainder of the paper. Hence, the complex three-vector (I29D can be 
re-written as 



± _ (- n afeM n 6a) = (- n afrM n fea)( n kJ " n ^ ' 
n ab\ n ba) \(- n ab\ n ta/ 



\ri \ "at I " 1 11 ba/ \ "ab\" l"6q/\"6ql 11 ab/ /'QfA 

V af>~ / +1 ±\ - 1/ _+i_+\i9 • \ 6U ) 



ab~ ( _ n ± |n ^ 

Next, we will need the following lemma 



Lemma 1. For all couples of coherent states (|ni),|ri2)) associated to the couple (111,112) of 

elements of S 2 , the following identity holds: 

(ni|ff|n 2 ) (n 2 |ni) = -(ni + n 2 - ini x n 2 ), (31) 

where the symbol 'x ' denotes the three-dimensional cross-product. 

Proof. The projector P n on the coherent state |n), can be written in terms of Pauli matrices 

P n = |n)(n| = -(1 + a-n). 

It obviously satisfies idempotency P n o P n = P n and unit trace trP n = 1. Hence, the left hand 
side of the above equation reads 

(ni|cr|n 2 ) (n 2 |ni) = - ((ni|<r|ni) + n 2 (ni|l|ni) + m 2 x (ni|<r|ni)) 

= ^(ni + n 2 - ini x n 2 ), (32) 



1 Here, we are supposing that (n a!) |X*J — nf, a ) 7^ 0, for all a < b = 1, 5. 
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where we have used (n|tr|n) = n and (n|n) = 1 in the last step. □ 

Using the above lemma, the complex three-vector (]30p can be split into real and imaginary 
parts, and the stationary equations (|28|) become the following twenty real vector equations 



E^ i n r^ =0, and n r^ =0, (33) 

b-.b^a 1 n ab' n 6a 6: b ^ a 1 U ab ' U ba 



for all o = l, .... 5. 



7 > 1 case. The variation for 7 > 1 proceeds in the same way but in this case the stationary 
equations are 

££V£ = ° (34) 

b^a 

with 



Vafe (m a6 |(X±)-lX±|m 6a ) ' 



plus variational equations for m which are vacuous at the critical points, discussed in the next 
section. 



5.1.2 Critical points 

As already remarked, the action (|15p is complex, and standard stationary phase does not apply. 
The key point is to supplement the stationarity condition with a maximisation condition on the 
real part of the action. The restriction of the stationary points to those maximising the real 
part of the action are called critical points. They contribute to the asymptotic formula while 
the other points are exponentially damped. The integrand of the four-simplex amplitude (jlOp is 
of the schematic form 

t „\a \ib 

1 ~ e e , 

where a < is the real part of the action and b is the imaginary part. When A — » 00, the only 
significant contribution to the integral comes from the configurations where e a = 1, i.e., when a 
hits its maximal value a = 0. 



7 < 1 case. More precisely, the real part of the action (I15p in the 7 < 1 case is given by 

Re 5 7< i = £ iab ln 5(1 " ^ • O +tib In 5(1 " "& • <J> (36) 

a<b 

where we have used the expression of the inner product between coherent states and all phases 
have been absorbed in the imaginary part of the action. The maximum Re S = of this 
expression is obtained when all n^ b and are anti-parallel, i.e., when for all a < b = 1, 5, 

Xtn ab = -X±n fea . (37) 

Thus, the dominating points in the integrand of the four-simplex amplitude (|l(jp are given by 
the points satisfying stationarity (I33p and maximising the real part of the action (I37p . All other 
points are exponentially suppressed. Evaluated on such points, the second stationarity equation 
becomes vacuous, and the critical points are those satisfying the ten vector equations (|3T1) . plus 
the five vector equations 

Kb n ab = (38) 

6: b^a 

for all a = 1, 5. We will refer to these two equations as orientation and closure respectively, 
for reasons that will become clear shortly. 
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7 > 1 case. The critical points for the 7 > 1 case are only slightly more complicated. The 
only significant contribution in the asymptotic limit comes when the following are satisfied 

— n ab = m ab 
m ba = n ba 

X+m ab = A"+ m ba 

X~m a6 = X~m ba . (39) 

for a < b. However, these conditions are equivalent to the set of equations (fHT|) . The action (JTHJ) 
at the critical points, when expressed as a function of 7 and the boundary data, gives the same 
function as (|15p but for the different range of 7. 

The integral coming from the resolution of the identity is no longer relevant as we assume 
the measures dm ab , dm ba are normalised to one. Using equations ([39]) in ([MJ) along with the fact 
that k ab = j~a b — j~ b , one can see that the stationary point conditions for 7 > 1 are equivalent to 
those for 7 < 1. 

5.2 Geometrical interpretation 

In this section, we show how geometrical structures emerge from the solutions to the critical 
point equations. 

5.2.1 Bivectors 

Let A 2 (M 4 ) be the space of bivectors. A pair of vectors N, M € M 4 determines a simple bivector 
NAM, and a general bivector is a sum of simple bivectors. A simple bivector B can be considered 
as an antisymmetric tensor 

nam = n®m-m®n, 

and hence an element of the Lie algebra so (4), using the Euclidean metric. This gives the action 
on a vector Q as 

(N A M)Q = (M • Q)N - (N ■ Q)M. (40) 
The norm \B\ of a simple bivector B is defined so that 

|JVAM| = |iV||M|sin0, (41) 

with 9 the angle between the vectors. Consequently, in terms of the components B IJ (for 
I , J = 0, ... ,3) of the antisymmetric tensor B this is 

\B\ 2 = -B IJ B U , 
2 

the indices being raised or lowered using the standard Euclidean metric 5. 

Let * be the Hodge operator acting on A 2 (IR 4 ); (*B) IJ = \e IJ KL B KL , where ejjKL is the 
totally antisymmetric four-dimensional tensor (g 0123 = l). With this definition, the Hodge 
operator satisfies ** = 1 and preserves the norm of a bivector. 

In the following it will be useful to split the bivectors into two three- vectors. Consider the 
eigenspace decomposition of the Hodge operator A 2 (M ) = A 2 _(M 4 ) A^_(]R 4 ) decomposing any 
bivector B into self-dual b + and anti-self-dual components, = ±b . We will extensively 
use the vector space isomorphism A^M 4 ) = R 3 explicitly given by 

6±i = ± S 0i 1 i B jk^ ( 42) 
2 J 
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with i, j, k = 1, 2, 3. Using the standard Euclidean inner product on M 3 , this gives 

|b+| 2 + 1 1» | 2 = 2\B\ 2 . 

This means that for a simple bivector B, |b + | = 1 1» | = \B\. 

If B and C are bivectors, then they form a Lie algebra with [-B,C] = BC — CB, using 
the action (|4(jp on vectors. The above isomorphism can be extended to an isomorphism of Lie 
algebras 

so(4) A 2 (M 4 ) = M 3 © M 3 = su(2) © su(2), 

using the map 

(b~,b+) i— ► (b~ • L , b + • L) 

and the notation b • L = b l L{ with Li the generators of su(2). A calculation shows that this is 
an isomorphism providing Lj = ^aj in the fundamental representation, so that 

[Li,Lj] = —€ij k Lk. 



5.2.2 The geometrical 4-simplex 

A geometrical four-simplex a is determined by the position of its vertices in M 4 . It is said to 
be non-degenerate if the five vertices do not lie in a hyperplane. In this case one also refers to 
the closed set formed by the convex hull of the vertices as the geometric 4-simplex; these can be 
glued together to form Regge manifolds. 

Each oriented triangle A of a geometric 4-simplex a determines a bivector B&{a) by wedging 
together two of its edge vectors. However the constructions presented in this paper naturally 
label the triangles of a 4-simplex by the two tetrahedra which intersect at the triangle. Suppose 
N a (o~) is the outward unit normal vector to tetrahedron a. Then N a (o~) AN b (o~) defines a bivector 
which is in the plane orthogonal to the triangle where tetrahedra a and b intersect. Therefore 
*(N a (a) AN b (a)) lies in the plane of the triangle; normalising it correctly then gives the definition 
of the bivector for the triangle labelled by the ordered pair ab. If k ab is the area of triangle ab, 
then this bivector is 

r f \ u N a (a)AN b (a) 

BM = k « b * \N a (a)AN b (a)\ - (43) 

This definition has the property that B ab (a) = —B ba (a). Also, the bivectors are unchanged 
if a is replaced by Ga, where G is an inversion x — > —x of M 4 , or a parallel translation. Apart 
from these symmetries, the bivectors determine the 4-simplex geometry uniquely. 

To state this theorem precisely, a general definition is required. An arbitrary set of ten 
bivectors B ab is said to be a bivector geometry if the following constraints are satisfied 

• Orientation. A bivector registers the orientation of the triangle to which it is associated, 
i.e., 

B a b = —Bba- 

• Closure. Each tetrahedron building up the four-simplex must close, 

Vo, Yl Bab = °- 

b: b^a 

• Diagonal simplicity. Each bivector must define a geometric plane, i.e., is a simple bivector 

B ab AB ab = O \b+\ = \b- b \. 
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• Cross simplicity. Each couple of bivectors belonging to the same tetrahedron must define 
a three-dimensional hyperplane spanned by the two planes associated to each bivector: 

Va, 3N a € M 4 , such that N aI B ! J = 0, V6 / a. 

• Tetrahedron. For three triangles meeting at a vertex e of the a-th tetrahedron (bcde), 

tr (B ab [B ac ,B ad }) /0, 
where we have used the isomorphism A(M 4 ) = so (4). 

• Non-degeneracy. The assignment of bivectors is non-degenerate. This means that for six 
triangles sharing a common vertex, the six bivectors are linearly independent. 

The following theorem is proved in [5], and more explicit detail is given on the tetrahedron 
constraint in [26j. One difference is that the constraints have been adapted to label the triangles 
by the missing vertices, which means that the bivectors depend on the orientation of the 4- 
simplex. 

The cross simplicity constraint stated here is as recently reformulated in pXJ. Note that the 
cross-simplicity constraint of [5] requires only that two bivectors B a b, B ac respect B ab A B ac = 0. 
This implies for simple bivectors sharing an index that they either belong to the same three- 
dimensional hyperplane or share a common direction. However in the case of the common 
direction the tetrahedron condition is not satisfied, as shown in pJJ. Hence the formulation given 
here is equivalent. 

Theorem 2. (Bivector geometry) The bivectors B a b{a) of a geometric ^-simplex a satisfy the 
bivector geometry constraints. Conversely, given a set of bivectors B ab on an abstract ^-simplex 
satisfying these constraints, there is a parameter /i = ±1 and geometric ^-simplex a, unique up 
to translation and inversion, such that B a i,(a) = ^B a ]y. 

We now show how the bivector geometry theorem relates to our critical point equations. 
We start by discussing the general idea. Consider an arbitrary bivector B = (b~,b + ). This 
bivector satisfies the diagonal simplicity constraint if and only if |b~| = |b + |. This is equivalent 
|13] to stating that there exists a Spin(4) rotation G = (X~,X + ) such that b^ = X^v, with 
v a fixed, arbitrary vector in M 3 . If G = (X, X) belongs to the diagonal SU(2) subgroup, the 
left and right handed components of B are equal as vectors b = b + and it is easy to see 
that a such bivector satisfies NiB IJ = 0, with J\f = (1,0,0,0) the S 3 'north pole'. In other 
words, B lives in a hyperplane orthogonal to J\f. In the general case, if G does not belong to 
the diagonal SU(2) subgroup, one can show that the bivector B = (X~ , X + )(v, v) satisfies the 
equation NjB IJ = 0, with N = GN ', and thus lies in a hyperplane orthogonal to the rotated 
vector N . To summarise, a bivector B satisfying the diagonal simplicity constraint can always 
be written B = (X~v, X + v) with the spin group element (X~,X + ) encoding the information 
on the normal defining the hyperplane in which the corresponding bivector lives. This idea will 
prove instrumental in the geometrical interpretation of the critical points. To start with, we 
introduce the following lemma. 

Lemma 2. Let a be a geometric four-simplex in R 4 with the areas of its triangles being half- 
integers. Then there exists a set of boundary data B, together with a set X^ satisfying closure 
(]38p and orientation A3TH . associated to a. These are detemined uniquely up to the action of 
SU(2) at each tetrahedron. 

Proof. Consider a geometrical four-simplex a. It is determined by five tetrahedra, each 
spanning a hyperplane to which we associate an outward normal N a (a) in S 3 , four out of which 
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are linearly independent, and ten triangles from which we construct the ten simple bivectors 

B ab (*) = (h- b (a),b+(°))- 

We introduce the diffeomorphism ( : S 3 -> SU(2); N i-> ((N) = N°l + iN^cXj between 

the space of unit normals S 3 C M 4 , and the unitary group SU(2). Prom this map, we can 

relate four-dimensional rotations G to the action of the spin group Spin(4) on SU(2): Q(GN) = 

X_C(N)X^. Using this one can see that given the five normals N a (a) one can pick five Spin(4) 

group elements G a = (X~,X£) via the transitive action of Spin(4) on S 3 : 

C(N a (a)) = X-C(AOPC)- 1 = X-(X+)-\ 

These are not unique, different choices being related by the diagonal SU(2) subgroup which 
leaves Af invariant. This acts separately for each tetrahedron and is the action of SU(2) in the 
statement of the lemma. 

Now pick one such set of values {X^ 1 }. From the bivectors B ab (o~) and the group element 
G a , construct the bivectors 

&c+) := ((X-r\{X+r 1 )(b- b (a),b+^))- 

From the discussion above, it is easy to see that this bivector lies in a hyperplane orthogonal to 
G~ 1 N a = Af. Accordingly, it necessarily has equal left and right components c~ b = c^ b := k ab n ab . 
(This is true for any such choice of {X^ 1 }.) The closure and the orientation of the bivectors 
B a b{°~) respectively imply, by inverting the above equation, the critical point equations ([37|) □. 

Note that the condition that the areas are half-integer is only required to meet the definition 
of boundary data. Thus the lemma would still be true if this condition were removed and the 
definition of boundary data did not require the areas to be half-integer. 

We can now state the converse, which is the key to interpreting the critical point equations. 

Theorem 3. (Reconstruction theorem) Given a fixed, non- degenerate boundary data B 
satisfying closure, J2 bj L a k ab n ab = 0, a set X^ 6 SU(2) ; a = 1, ...,5, solving the orientation 
equations X^n ab = —X^n^a is one of two geometric types: 

• The solutions determine a non-degenerate geometric ^-simplex a in M 4 , defined up to 
translation and inversion, and a parameter fj, = ±1, or 

• The solutions determine a single vector geometry. 

In the first case, the bivectors B a b(a) of the ^-simplex a satisfy 

B a b{<y) = nkab{X~ , X+)(n ab ,n ab ), 

for all a^b. 

Proof. The proof is based on the bivector geometry theorem. Consider the set of ten bivectors 
B ab constructed as follows: 

Bab = (b^b+J := k ab (X~ , X+)(n ab , n ab ). 

These satisfy the orientation, closure and diagonal simplicity bivector geometry conditions. The 
cross-simplicity condition follows from the fact that the bivector (n ab ,n ab ) is orthogonal to the 
north pole Af = (1, 0, 0, 0), and so the four bivectors B ab with fixed a are orthogonal to the vector 
iV a obtained by the action of (X~ , X^ ) on Af. The tetrahedron condition follows from the fact 
that tr (B ab [B ac , B ac i]) is proportional (by a non-zero constant) to plus or minus the squared 
volume of the geometric tetrahedron determined by the boundary data, which is non-degenerate 
by hypothesis. 

We have therefore constructed bivectors satisfying all of the bivector geometry constraints 
except for the non-degeneracy condition. To investigate this final constraint we need the following 
lemma. 



19 



Lemma 3. Given a boundary data setB satisfying (|37p . the SU(2) group elements satisfying 
(I38p are such that the corresponding constructed five unit vectors N a £ S 3 fall into one and only 
one of the following two classes: 

1. At least three out of the five N a are linearly independent 

2. The N a are all proportional, i.e., N a = ±N for some fixed N. 

In the first case the constructed bivectors satisfy the non- degeneracy condition of theorem^ 
To establish the first part of the lemma we show 

• If three of the normals {N a , N b , N c } are proportional to N , then consider a fourth normal 
N<i. By the gluing condition the bivectors B ea , B e )j, B ec lie in the hyperplane orthogonal 
to N and so by closure so does B e d- Thus the tetrahedron lies in this hyperplane and so 

is also proportional to N. The same holds for N e . 

• In all the other possible cases, three of the normals {N a , N b , N c } must be pairwise linearly 
independent. Then one can conclude that all three together are a linearly independent set. 
This follows from the fact that the intersection of hyperplanes a and b and the intersection 
of hyperplanes a and c are non-parallel planes in hyperplane a (due to the fact that they 
are faces of a non-degenerate tetrahedron). 

Now, given three linearly independent normals to tetrahedra {N a , N b , N c } we will construct 
six linearly independent bivectors from the tetrahedral bivectors. We take x^x 2 ^ to generate 
the plane H iV^ and from these two vectors we complete with xa and with xb , we 
do this in such a way that xb also completes n , which can always be done. Then we 
construct the five linearly independent bivectors 

{ x lb A x lb> x lb Ax A, xl b Ax B , x 2 ab Ax A , x 2 ab Ax B }. 

Each of these lie in the hyperplane of a tetrahedron, and thus is a linear combination of the 
terahedral bivectors. Together with xa A xb they span the whole space of bivectors. 

Now we take the intersection Njj- n generated by x bc , x\ c and complete with xc- We 
have that xc = ctix l ab + oi2X 2 ab + a^xs + a^XA since these span the whole 4 dimensional space. 
Furthermore «4 7^ since otherwise and Njj- would be spanned by the same vectors and 
N c || N b . Furthermore x\ c = (3ix^ b + ^2x\ b + 03Xb can be chosen such that fa is non zero, 
since n N b ^ N b PI , due to the non-degeneracy of the boundary data. Then the new 
bivector x\ c A xc contains the nonzero term a^P^XA A xb along with terms in the set of five 
bivectors above. Therefore a complete set of six linearly independent bivectors is generated by 
the tetrahedral ones, and Lemma 2 is proved □ 

The above lemma implies that the set of constructed bivectors B ab are one of only two 
types. In one case they satisfy all bivector geometry constraints, including the non-degeneracy 
condition. Then using theorem [21 we can conclude that there exists a parameter \x = ±1 
and a geometric 4-simplex a, up to translation and inversion, with bivectors B ab (o~) given by 
B a b(o~) = fJ-B ab , determined by the corresponding solutions to the critical point equations. 

If the constructed bivectors do not satisfy the non-degeneracy condition, they all lie in the 
same hyperplane and hence determine a vector geometry which is not a non-degenerate geometric 
4-simplex. This concludes the proof of the reconstruction theorem □ 

5.3 Classification and uniqueness of the solutions 

An important input in the asymptotic formula is the classification of the solutions to the critical 
point equations. 
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5.3.1 Symmetries of the critical points 

To start with we need to consider their symmetries. These include the symmetries of the action 
discussed in 13.41 Given a set of five SU(2) elements {U a } solving (|37|) we immediately see that 
{e a YU a } with Y G SU(2) and e a = ±1 is also a solution. The sign ambiguity is due to the fact 
that two SU(2) elements differing by a sign act in the same fashion on three dimensional vectors. 
If two solutions are related by the symmetries we write {U^} ~ {U^}. If they are not related by 
symmetries we write {U^} oo {U%} and call them distinct. 

Given two solutions to the critical point equations we have two ways of forming Spin(4) 
elements (X£,X~) out of them, that is X^ = f7„, X^~ = U%. Thus we also have 

P: (X-,X+)^(X+,X-) 

as a symmetry. This is not a symmetry of the action but a symmetry of the solutions of the 
critical point equations. 

The induced action of Spin(4) and P on the normal vectors N a = £ _1 (X~(X£)~ are the 
action of SO (4) and the parity transformation 

P G 0(4), P : («o,v) h-> («o,-v) where (vq,v) G (RAT) eTV" 1 

respectively, acting on N a G R 4 . Together, these generate the group 0(4), the full set of linear 
isometries of four-dimensional Euclidean space. 

We can now translate the result of lemma [3] into a statement on the group elements of the 
solutions. Note that using the isomorphism between S 3 and SU(2) we can turn the results 
for the normals in the degenerate sector, N a = ±iV, into X~ (X+) _1 = e a Y. Thus we have 
X~ ~ X£. The converse clearly also holds. Since we showed in lemma [3] that the only other 
case is full non-degeneracy we also have that if the two sets are distinct these group elements 
define a non-degenerate 4-simplex up to inversion. 

Bivectors register the S0(3) x S0(3) quotient of the Spin(4) symmetry. The parity trans- 
formation P flips the selfdual and antiselfdual sectors when acting on bivectors: P(b _ ,b + ) = 
(b + , b~). This can be seen from the definition of the sectors (|42[) : 

P : ±B 0i + ^e l jk B jk ^ ^B 0i + ^e l jk B jk . 

The second term has two spatial indices and so does not change sign. 

As bivectors do not register the full SU(2) symmetry there is a residual set of symmetries 
not acting on them given by the signs e„ . Part of this symmetry is however registered by the 
normals used in the reconstruction theorem. These transform as X~ (X+)~ l — ► e^e~ X~ (X^)^ 1 , 
and thus N a — ► e^e a ■N a . 

In the case of a non-degenerate 4-simplex, the reconstruction theorem gives a pair of 4- 
simplexes a, a' related by inversion. The symmetries can be used to turn the tetrahedron 
normals into ones which are all outward pointing for one of these and all inward pointing for 
the other. At such a critical point then a and a' are distinguished, and the one which has all 
normals outward pointing will be used as the canonical 4-simplex associated to the critical point. 
This 4-simplex is denoted a in the following. 

The symmetries do not affect the 3- vectors n a {,. However they do affect their relation to 
the 4-simplex a. For a given 3-simplex a the inverse of the SO (4) transformation determined 
by X~,X£ maps this face of a to a tetrahedron t' a in M 3 (the hyperplane orthogonal to J\f). 
The n a b are either all outward-pointing normals for the faces of t' a , or they are all inward- 
pointing normals. The orientation condition implies that if they are outward (inward)-pointing 
in one tetrahedron, they are outward (inward)-pointing, respectively, in all the others. In fact, 
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a calculation shows that the parameter \x = ±1 of theorem [2] for the 4-simplex takes the value 
+ 1 for outward-pointing tetrahedron normals and —1 for inward-pointing. 

The action of P flips [i — > — fj,, whereas the action of Spin(4) leaves it unchanged. This 
follows from the fact that PB ab (o~) = —B a b(Pa), as can be seen from the presence of the Hodge 
star in equation (|43p . Thus ji measures the orientation of the 4-simplex; [i = 1 if the orientation 
agrees with the orientation of its boundary specified in section [3.3.21 and [i = — 1 if it has the 
opposite orientation. 

5.3.2 Classification 

We can now prove the following lemma: 

Lemma 4. Let B denote the set of fixed boundary data {Ta. ab , & a &}a^6=l,...,5j with n in S 2 , and 
k £ N/2 satisfying closure 

Kb n a6 = 0, (44) 

b: b^a 

for all a = 1, ...,5. Consider a set of five SU(2) group elements U a solving the equations 

U a n ab = -U b n ba . (45) 
There are at most two distinct sets ofU a . 

Proof. Assume there are three distinct solutions to (j4"5|) . say {U^}, {U% }, {U® }• Then we 
can form non-degenerate solutions to the critical equations, e.g. (X~,X+) = (U^,U^) and 
(X~, X£) = (U%, U®). These then give two non degenerate geometric 4-simplices with the same 
boundary data up to inversion. The boundary data however gives the full set of tetrahedra 
forming the boundary of the geometric 4-simplex. However any two n-simplices with congruent 
(n — l)-faces are in fact congruent themselves (as discussed in section [3]). Therefore the two 
4-simplices must be the same up to rigid motion. In particular they are related by P and 
SO(3) x SO(3) up to inversion. Thus we have either that {U^} ~ {U%}, or, if the rigid motion 
mapping them to each other involves a parity, that {U^} ~ {U®} and {U®} ~ {U 2 }. Both are 
in contradiction to our assumption that {U*},{U%},{U%} are distinct. It follows that we can 
never have three or more distinct solutions to (|45jl .D 

From this it follows that there are three types of B, distinguished by how many solutions 
{U a } to flM!) they allow: 

1. The set B is such that there are no solutions to (fl5l) . This means the tetrahedra do not 
glue at all. 

2. There is only one distinct solution to (f4"5j) . thus there is only one set of solutions to the 
critical point equations. This is given by (X~,X£) = (U a , U a ). 

3. There are two distinct solutions to (|4"5|) and we have, up to the symmetries of the action, 
four solutions to the critical point equations: 

(X~,X+) e {([/-, U+), (U+, U~), (u+, u+), (u~,u-)} 

These are in one to one correspondence to non degenerate geometric 4-simplices up to rigid 
motion. 
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5.4 The action at critical points 

In this section the action at a critical point is expressed in terms of the underlying geometry. 
The critical points satisfy closure and the conditions 

X±{n ba ) = -n a6 , (46) 

for all a ytz b. The lift of this equation to the coherent states involves a phase 

X±\n ba ) = e i ^\-n ab ). (47) 

The phases 4>^ b and <fi~ b are defined modulo 2tt. However it is the sum and difference of the 
phases that enters into the geometrical expressions. 

Evaluating the action S = iA, either (|15p or (|18p . on the critical points leads to 

A = J2 k ab (4>+ + 0- ) + 7 k ab (0+ - 0" ), (48) 

a<b 

for all 7/I, where we have used the corresponding relation between and A;. In this formula, 
the phases enter via the variables 4>~£ b + <^>~ fe and (j)^ b — (f)~ b which are both defined mod 4-7T, but 
with one remaining indeterminacy given by the simultaneous replacement 

^b+ ( t>ab^4> + a b + Kb + ^ 

Taking into account the fact that the k ab are half-integers and that (1 + ^)k ab are integers, one 
can see explicitly that A is well-defined mod 2ir. The indeterminacies may be fixed by choosing 
suitable coordinate ranges; we will use 

-2tt<0+ + ^<2tt, 

Further analysis of the action is done by separating two special cases. 
5.4.1 Non-degenerate 4-simplex, X + and X~ distinct 

In this case, the critical point is given by the sets {^q^} and {^^} which are assumed to 
be distinct from each other, i.e., not related by a symmetry. According to theorem [3l these 
correspond to a pair of geometrical 4-simplexes a and a 1 , related by inversion. The results of 
this section will show that the phase difference term in the action (I48p evaluates to the Regge 
action of the geometric 4-simplex. Further, in the case of a Regge state it will be shown that 
the phase sum term in the action vanishes. 

According to the discussion in 15.3.11 there is a discrete symmetry of the action which relates 
any critical point to one in which N a € M 4 is the outward pointing normal of the 4-simplex a for 
each tetrahedron a. Therefore it is enough to evaluate the action on these critical points. 

The Regge action for a flat Euclidean four-simplex is a boundary term that involves the 
dihedral angles < G < tt between the normals of the boundary tetrahedra [6|. This angle is 
calculated through the scalar product between the outward normals to the tetrahedra 

cos Q ab = N a - N b 

= \ lr{£{N a )C{N b y^ 

= \ ^(X- b (X+)- 1 ), (49) 



23 



where we have used the fact that we are considering only the critical points corresponding to 
outward pointing constructed normals, that is, N a (o~) = N a for all a. The key point is that the 
above angle can be related to the rotation angle associated to the rotation stabilising the vector 
n ab , since the coupling of two of the equations (1l6j) leads to the relation 



KbiKbr^ab = n ab . (50) 

Similarly, (|4"7|) leads to 

X o - fe (X+)- 1 |n o6 ) = e^-^|n a6 ) (51) 

This shows that the SU(2) matrix X^X^,) -1 generates rotations about the n ab axis, and has 
eigenvalues e^af, - *^) and its complex conjugate. It can be written 

X d>( X tb)~ l = ex P - <Kb)° ■ n ab = exp 2(</>+ - 0- )L • n ab (52) 

using L = as before. 

Equation (jl9j) then shows that 

\Kb - <Kb\ = Q ab- 

This argument does not resolve the sign of the left-hand side, and a more delicate argument is 
required. The core of our argument is that this sign depends only on the relative orientation of 
the bivectors *B ab and the geometric 4-simplex, in fact: 

Lemma 5. The relative sign between the angle 4'ab~ ( l ) ab an d ^ e dihedral angle Q ab is determined 
by the parameter fj, = ±1 controlling the relative sign between the geometrical and constructed 
bivectors of theorem 0' 

<^ab " <Kb = V @ ab- 

Proof. The proof is based on the explicit construction of the dihedral rotation associated 
to the dihedral angle. Recall that a simple bivector B ab admits two stabilising Spin(4) sub- 
groups. This is due to the isomorphism A 2 (M 4 ) = spin(4) and to the fact that spin(4) admits 
a two-dimensional Cartan subalgebra. This is reflected in the fact that B a b, regarded as a Lie 
algebra element, obviously commutes with itself but also with its Hodge dual *-B a t>- Hence, the 
Spin(4) subgroups G ab and G^ b generated respectively by the Lie algebra elements -B a fe/|-E>ab| 
and *B a t,/\B a t l \ stabilise B ab . Geometrically, the subgroups G ab and G^ b correspond to simple 
rotations acting in the plane defined by B ab and in its orthogonal complement respectively. 

Consider a geometrical 4-simplex a and a particular bivector B ab {a\ The corresponding 
isotropy group G^ b contains an element of particular interest; the dihedral rotation D ab £ SO (4) 
mapping the normal N a to the normal iVj and stabilising the orthogonal plane 

3ab = exp i Qab TNWk\) ■ (53) 

This is clearly a simple rotation in the plane spanned by iV a and N b , so we just need to check 
that the rotation has the right direction and scale factor. To do this we expand the exponential 
above to first order and apply it to iV a obtaining 

-n m n Nb A N ° \ at m j- ®ab(N b -N a cos Q ab ) 

1 + Qafe | AT ~ AT , • N a = N a + ; . , (54) 

\N b AN a \J |sm6 ab | 

where we have used (|4ip and (|40p . For small positive @ ab , clearly the expression above reads 

D ab ■ N a <* N b , 
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which means (J53J) is correct. 

Using the definition (I43h , theorem [3] and the isomorphisms of section 15.2. 1\ the bivector is 
expressed in terms of the boundary data as 

|jy° A jy b | = ~K X a ' x +)(n ab , -n ab ) = (-fiX~n ab , fj,X+n ab ). 

Introduce the exponential 

dab = ±exp nQ ab n ab • L e SU(2). 
A lift of the dihedral rotation to SU(2) x SU(2) is 

D ab = (X->d ab ,X+>d~ b 1 ) (55) 

where the symbol '>' denotes the adjoint action of SU(2) on itself, X > d = XdX^ 1 . Either 
choice of sign in d ab gives a lift, reflecting the 2-1 ambiguity in the spin group. The choice of 
sign is arbitrary for the purposes of this proof. 
Finally, the equation D ab N a = N b implies that 

C(D ab N a ) = (X~ > exp/i9 ab n afc • L) X~(X+)~ l (X+ D> exp -/iG afe n a6 • L)" 1 

= X- (exp 2^9 a6 n a6 • L) (X+)" 1 = X b (X+)' 1 (56) 

Using ([52]) now gives 

exp 2fj,Q ab n ab • L = exp 2(<^ - <p~ b )n ab ■ L 

and hence n@ a b = <Pab ~ ^ab- This concludes the proof of lemma □ . 

In rest of this section it is shown how the action (I48D simplifies for the case of a Regge state. 
In ([55]) a key role is played by the SU(2) element 

dab = K a b exp flQabTO-ab " L, (57) 

with K a b = ±1. By (|52p we have 

V,/, 

Let us now define 

lab^ab aI1Q n ab : = a ab^ab 

We can check in fact that k^ b = h~ b =: h ab by calculating (h~ b )~ 1 h^ b = (h^ b )~ l h~ b = 1. 

The key point now is that the rotation (X~ b , X^ b ), that maps M ± into Af ± , can then be 
uniquely decomposed into two simple rotations: 

K> Kb) = ( d ab, d-J) (h ab , h ab ), (60) 

The element (h ab ,h ab ) belongs to the diagonal SU(2) subgroup of Spin(4) leaving M invariant 
and is thus a rotation in the hyperplane M = M 3 . This three-dimensional rotation glues the 
triangle corresponding to the bivector (n ba ,n ba ) onto the triangle associated to the bivector 
(n a b,n ab ). Note that by (|47p . we have: 

d~ b h ab \n ba ) = X+\n ba ) = e 1 ^ J \ n ab } 

dabhab\n ba ) = X~ b \n ba ) = e 1 ^ J\n ab ) (61) 



dl = X- b (Kb) 1 (58) 



Kb ■■= dabX+ and h~ b := (59) 



Now consider the following diagram: 



t a {X - X »\ T a (62) 



{g a b,g a b) 



D, 



a b 



h — n 
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where t a € M 3 C M 4 are the tetrahedra at the north pole (i.e. in the hyperplane orthogonal to 
TV), r a G M 4 are the actual geometrical ones in the 4-simplex a and the g a b £ SU(2), defined in 
Q, are the maps that effectively glue the triangles 4> a (A a b) and </>&(A a 5), and leave the normals 
inside t a and tb anti-parallel. Note that by the reconstruction theorem, the maps in the diagram 
commutes when acting on both the triangles (f) a (A a b) and on the internal normals n a ;,. Since 
furthermore all maps in the diagram are orientation preserving, the composition maps coincide 
on a three dimensional subspace and are orientation preserving, which means the SO (4) action 
of the maps in the diagram commutes. This means the diagram of maps as elements of Spin(4) is 
commutative up to sign. Therefore it is possible to now determine the sign K a b in the definition 
of (1^ as the one that makes this diagram of maps commutative in Spin(4). 
Hence we have that 

( X adab, X £d~b ) = i X b 9ab> X b 9ab) 

which gives 

9ab = { x b )~ X~d ah = (X b )~ x +d~ b , 

and which also tells us that h a b = g~ b . If we substitute this fact in (I6ip and use the action of 
g~ b on the coherent state \rib a ), then we obtain 

d- b l J\n ab ) = e^J|n a6 ) 

dabJ\n ab ) = e^J|n a6 ) (63) 
Taking the eigenvalues, and writing n a b = (— l)" ab with v a b = or 1 gives 

Kb = \^®ab + KV ab 
1 ^ 

Vab = ~2^ab + KV ah 

mod 2tt. Adding these gives 

Kb + <t>ab = 27TU ab 

as an angle between and 4-7T. 

The signs K ab cannot be always equal to one, since the symmetries g a b — > e a ^b9ab, for e a = ±1, 
change them by 

K ab l—> £a,£b K ab- 

Exactly the same transformation also occurs as a result of applying the symmetries (X~ , ) *— * 
(e a X~,e a X+). ^ 

However it is possible to show that these symmetries can be used to transform either the g 
or the X so that n a b = 1 for all a and b. 

Lemma 6. The signs defined in (|57p obey K a b = e a eb for some e a = ±1. 

Proof. The geometrical interpretation of the dihedral rotation D a b € SO (4) is that it restricts 
to the Levi-Civita connection u a b on the boundary da. Let r a be the a-th tetrahedral face of a. 
Then the dihedral rotation maps the tangent spaces 

D ab : Tr a ^TTb. 

Therefore the lift D a b € Spin(4) represents the spin connection in some gauge; indeed according 
to (f62|) it is gauge-equivalent to the g a b, which does indeed represent the spin connection on the 
boundary. 
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According to the discussion in section 13.3.31 the holonomy of the spin connection around an 
edge in the boundary 3-manifold da is uniquely determined; it is a rotation by the deficit angle 
of the Regge geometry, modulo An. This can be characterised as the rotation which deforms 
continuously to the identity rotation (and not a 2tt rotation) when the geometry around the 
edge is deformed to the flat geometry with zero deficit angle. Alternatively, one can view this 
as a statement about the contraction of the loop to a point in a smoothing of the geometry. 

The Levi-Civita connection has a canonical lift, given by 

d ah = exp jlQabTO-ab ■ L = Kabdab, (64) 

which is used to give the canonical lift of the dihedral rotation, 

D ab = {X- > d ab , X+ > d£). 
The holonomy around an edge of the boundary which lies in the 

D C aDb c Dab- 

This is the correct holonomy because it lifts the Levi-Civita connection and reduces to the 
identity element of Spin(4) if the dihedral angles 6 a ;,, 6b c , ca are deformed continuously to 
zero. Therefore this holonomy is equal to 

D ca DbcDab- 

Thus around each edge one has the cocyle condition 

l^ca^bc^ab — 1> 

which, since H 1 ^ 3 , Z2) = 0, implies that = e a et, for some values of e a = ±1. □ 
Finally, putting these results in the action (|48p results in 

A = fl^^kabQab. (66) 

a<b 

5.4.2 Non-degenerate 4-simplex, X + and X~ related by symmetry 

In this case {A"+} ~ {X~}, which implies that 0i = (f)~ b . Hence the action is 

a<b 

According to section 15.3.21 for a non-degenerate 4-simplex there exist two distinct solutions to 
the critical equations, {U^} and {U~}. The case here arises when {X+} ~ {X~} ~ for 
example. Putting X' a + = and X'~ = U~, then the reconstruction theorem applies to the 
data {X' a + }, {X' a ~}, as in the previous section. Then the results of the previous section applied 
to {X'+}, {X'-} show that 

Kb = 2^®ab + ™ ab . 

The parameters [i and Q a b are those determined by {X' a + }, {X' a ~}. 
Thus the action is 

A = k ab@ab, (67) 

a<b 

the sum over the f a & terms again vanishing by lemma El 



(65) 

three tetrahedra a, b and c is 
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5.5 Proof of theorem Q] 



The results of the previous sections can now be assembled into the proof of the main theorem. 
The idea of the proof is to apply the method of extended stationary phase to the integral formula 
(]14p for 7 < 1 and (I19p for 7 > 1. However before this is done, it is necessary to account for the 
symmetries of the formula. 

According to the results of section 15.3.21 given fixed boundary data, the number of distinct 
classes of critical points not related by the symmetries of the action (section |3.4.1|) is either none, 
one or four solutions. 

For any critical point, the discrete symmetries give 2 10 critical points which give the same 
contribution to the asymptotic formula. However four of these discrete symmetries (with 
independent of a) are in the class of continuous symmetries. Thus in each symmetry class there 
are 2 8 connected manifolds of critical points, each of which has the free action of the Spin(4) 
continuous symmetries. In fact in the case of a non-degenerate 4-simplex, a group of order 2 4 
was used to make the normals N a outward-pointing, and another group of order 2 4 was used 
to set K a b = 1. However it is actually not necessary to keep track of this: the contribution of 
each of these connected components is the same in each symmetry class and can be evaluated 
by looking at any one of them. 

In order to apply the method of extended stationary phase we need to ensure that the sta- 
tionary points are isolated. This is accomplished by the change of variables X^ a = (X^ 1 )~ 1 X^ . 
After this X^ drops out of the action, and we are left with an SU(2) x SU(2) integration which 
is redundant. The 'gauge fixed' integral formulas then have isolated critical points related only 
by the discrete symmetries and can now be evaluated using extended stationary phase. 

Consider first the case of four solutions. This is exactly if the boundary data corresponds 
to the boundary of a non-degenerate geometric 4 simplex. For this boundary data, consider the 
symmetry class of critical points with distinct X + and X~ . In section [5.4.11 it was shown that 
the action at these points evaluates to (|66|) . From section [5.3.21 we know that exactly two such 
solutions occur and that they are related by parity. From section 15.3.11 we know that parity only 
changes the sign of [i in (I66p . We can express the asymptotic result in terms of the solution 
{A~, X + } = {U~,U + } which has /i = —1. The phase of this asymptotic contribution to the 
amplitude is then given by 



with Q a b the dihedral angles of geometric 4-simplex corresponding to the boundary. The con- 
tribution from {U + , U~} has phase 



The two additional terms we obtain in section f5.3.2l correspond to the solutions {U + ,U + } 
and {U~,U~}. In section 15,4.21 we evaluated the action at these points as (|6T|) . Consider 
the {U + ,U + } solution. The [i occurring in the evaluation of its action is the same as in the 
evaluation of {U + , U~}. Therefore we have \x = 1 for this point and /i = —1 for {U~ , U~}. The 
phases of their contribution to the asymptotics are therefore 




a<b 




a<b 




a<b 



and 




a<b 
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respectively. 

As the final element of the asymptotic contributions of each of the solutions we need to 
calculate their Hessian matrices which are given by 



= I SOT 1 <68) 



for 7 < 1. This is now a 24 by 24 matrix as we have factored out the integration with re- 
spect to X 5 and c and d range from 1 to 4 only. We write H({X + , X~}) for the Hessian 
evaluated at the critical point {X + , X~}. The Hessian at the critical points is given explic- 
itly in terms of geometric data corresponding to these points in section 15.61 The expression 
for 7 > 1 is only marginally more complicated and is given there as well. In both cases we 
have det(if({X + , X~})) = det(H({X" , X + })) (see next section). This allows us to combine 
the asymptotic contributions from the two solutions with distinct elements into one cosine. The 
factor a(xo) in the stationary phase formula is the normalising factor for the measure on 8 copies 
of SU(2) in coordinates which agree with the derivatives; each copy is a factor of (47r)~ 2 . 

As each of the critical solutions occurs 2 8 times we thus find that the complete asymptotic 
behavior of (|14|) is given by: 



h(i>\) ~ (-i) 



V A J (4tt) 16 I y/det H({U+,U~}) 



exp(iA ^2 a<b k ab e ab ) exp(-«A £ a<6 k ab Q ab ) 



^det H({U+,U+}) ^det H({U-,U-}) J ' 
This establishes the non degenerate boundary data with 7 < 1 part of theorem [T] with 

N ^ 2 = ^(detH({X^,X^}))-h 

and ei, 62 = ±. 

For 7 > 1 we have four additional variables, m ab , m ba per propagator in the action. Fur- 
thermore the resolution of the identity contributes a factor of (d-+ ) 2 = (2j~t + l) 2 . Therefore 
we get an overall additional factor of 

(-l)Or-i)*- ( MA2(7 + l)fca6 + l) y ^ ( _ 1)(7 -i) fca6((7 + 1)kab) 2 
\ 47rA / 

per propagator. Thus the numerical factors for 7 > 1 are 

N7 ^ = T^( detH tt xei > xt2 ^ 2 II(- 1 ) (7-1) ^((7 + i)^ 6 ) 2 - 

^ ' a<b 

For a set of boundary data which is not a non-degenerate 4-geometry but forms a vector 
geometry, we only have one contribution to the asymptotic formula up to symmetries. For this 
case we did not evaluate the phase of the single contribution. The stationary phase formula 
gives 

/2tt\ 12 2 8 1 

|««| ~ (-J wF ^mmm\ m 

for 7 < 1 and the same prefactors as in the first case for 7 > 1. 
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For a set of boundary data that is neither a non-degenerate 4-geometry nor a vector geometry, 
there are no stationary points of the action and the stationary phase formula states that f$ 
decreases faster than any power of A. 

Ufa) = o{\- N ) ViV (70) 

This concludes the proof of theorem [T] apart from the facts used about the Hessian which 
are established in the next section. □ 



5.6 The Hessian 

This section provides an explicit formula for the Hessian matrix which appears in the asymptotic 
formula. The symmetry of the Hessian under the interchange of + and — is then apparent. We 
will begin by discussing the 7 < 1 case. Here the Hessian matrix of the full action is a 24 x 24 
matrix defined as 

FS±± = {dx^dxi*) (71) 

As the action is a sum of left and right sector, the terms off diagonal with respect to the ± index 
are and H is a block diagonal matrix 

where H ++ ,H are the 12 x 12 Hessian matrices for the left and right sectors respectively. 

As the Hessian is block diagonal, its determinant can be written as the product of the 
determinants of the two submatrices 

det(il) = det(H-) det{H ++ ). (73) 

For clarity we just give the result for det(H ++ ), det(H ) is the same. From (|15p . using again 
(I27p for the variation of SU(2) elements we obtain for the elements diagonal in the vertex index 
c: 

H % ++ = dx SL, + = -\ £ - - ^ + ) ( 74 ) 

c c b^c 

Where was defined in (|29p . Now if we denote a set of critical points by {X ,X + }, and 
use Lemma [T] to evaluate the Hessian at the critical points, we have 



(75) 



which is equivalent to the Hessian for the asymptotics of a coherent tetrahedron that was derived 
in [p2]. For the off-diagonal terms we obtain: 



d 2 s 



dxTdxi; 



{X+} 



l2j+(sV-U* k n* + -n%T#) (76) 



For 7 > 1 there is an additional dependence on the coherent state vectors m a b. This means the 
Hessian, while still symmetric in the left and right sectors, is now no longer block diagonal. 
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6 Remarks and Conclusions 



We have studied the semi-classical limit of the four-simplex amplitude of the Euclidean EPRL 
model for both sectors of the Immirzi parameter. Using the method of coherent states, we 
have shown that non-trivial contributions come from boundary states which determine a Regge 
geometry. By contrast, most non-Regge geometries have suppressed asymptotics. 

In the case of a Regge state associated to a non-degenerate 4-simplex boundary geometry, 
the asymptotic formula contains the cosine of the Regge action 



with the areas of the triangles given by "fk a b, as calculated in the original EPRL paper. Note 
that it is also logically possible to identify the areas with the spins k a b, in which case the Regge 
action appears in the argument of the complex exponential terms. Interestingly, these terms 
scale with the same exponent. 

We have also studied the semi-classical limit of the model for more singular boundary data 
that we have labelled vector geometries. For this type of boundary, the asymptotic formula 
is not dominated by non-degenerate four-simplex geometries and does not contain the Regge 
action. For fully non-Regge boundaries, the amplitude is suppressed in the large spin limit. 

An interesting corollary of our work is the asymptotic analysis of the EPR model where the 
Immirzi parameter is set to zero. In this case, the left and right spins are equal = j~ b = j ab 
and our result leads to the following asymptotic formula for a Regge-like boundary state 



The completion of this work required an appropriate posing of the asymptotic problem. This 
has led us to unravel the crucial role of the phase information contained in coherent states. We 
have indeed shown that, while the phase of coherent states is unimportant in the description of 
a single coherent tetrahedron, it plays an important geometrical role in the description of more 
complicated triangulated three-manifolds. We have discovered that the phase of coherent states 
registers the gluing properties of three-dimensional tetrahedra. Hence, the framework underlying 
our semi-classical analysis has a flavor of the area/angle formulation of Regge calculus by Dittrich 
and Speziale {32]. 

Several extensions of this work require attention. Firstly, the asymptotics for the Lorentzian 
models [T5J EI] should be computed. Secondly, a numerical test of the asymptotic formula in 
comparison to direct calculation of the spin network vertex amplitudes is desirable. Finally, a 
study of the renormalisation of the model and its behaviour under Pachner moves would be an 
interesting avenue of investigation. 
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